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Impossible role of ^He impurities in solid ^He 

o ■ 

O ; 

CNefstratios Manousakis 



5-H 

2 



< 

cn 
(N 



o 

00 

> 

m 
oo 
m 

o 

o 



Department of Physics and MARTECH, Florida State University, Tallahassee, FL 32306-4350, USA and 
Department of Physics, University of Athens, Panepistimiopolis, Zografos, 157 84 Athens, Greece. 



PACS 67 . 80 . -s - Solid helium and related quantum crystals 

PACS 67.80.Mg - Defects, impurities and diffusion 

PACS 67 . 90 . +z - Other topics in quantum fluids and solids 

Abstract. - We use a quantum lattice gas model to describe essential aspects of the motion 
of *i?e atoms and of "^He impurities in solid "^He. This study suggests that ^He impurities 
bind to defects and promote '^He atoms to interstitial sites which can turn the bosonic quantum 
disordered crystal into a metastable supersolid. It is suggested that defects and interstitial atoms 
are produced during the solid *i/e nucleation process where the role of ^He impurities (in addition 
to the cooling rate) is known to be important even at very small (1 ppm) impurity concentration. 
It is also proposed that such defects can form a glass phase during the *-H"e solid growth by rapid 
cooling. 



I 

o 
o 



X 



Kim and Chan [1] (KC) using the torsional oscillator 
technique found a decrease in the resonant period of solid 
ii?e confined in porous vycor glass and in bulk solid he- 
lium below 200mK, indicating the possible onset of su- 
perfluidity in solid helium. The experimental results of 
KC have been independently confirmed by Rittner and 
5-eppy [2] (RR) using a different geometry to confine the 
solid. In addition, RR observed that the superfluidity of 
^olid helium can be significantly influenced or eliminated 
by annealing of the solid helium sample. Because of these 
ijistory-dependent results and the negative results in at- 
tempts to drive flow by pressure [3], the interpretation of 
t]he results of KC is subject to debate. 
■ Furthermore, it has been observed that the superfluid 
r,esponse and superfluid fraction are strongly dependent on 
the amount of isotopic ^He impurities which exist in the 
' naturally available helium. When, for example, the natu- 
rally occurring concentration of ^He impurities of about 
one part per million is reduced to less than one part per 
billion [4] the superfluid fraction is reduced from about 1% 
to approximately 0.03%. When the ^He impurity concen- 
tration is increased to about 0.1% the superfluid fraction 
vanishes [1]. Previously, Ho et al. [5] found an anoma- 
lous behavior of the acoustic attenuation and velocity in 
solid ^iJe below ^ 2{)QmK at a low concentration of ^He 
impurities. 

The possibility of superfluidity of solid ^He has been 
extensively discussed [6-9]. It is of fundamental value for 
our understanding of this and related phenomena because 



of the implication that there can be coexistence between 
spatial and momentum space order [10]. 

In this paper we study the role of very low density of 
impurities in a bosonic hard core solid. While the theo- 
retical studies of the role of impurities in solid helium has 
a long history [6, 11], the present study, based on an anal- 
ogy with models of doped quantum antiferromagnets and 
by using rather recently developed techniques for such sys- 
tems, sheds some "new" light on the problem. Our model, 
that describes the impurity motion in an otherwise ideal 
quantum bosonic crystal, maps to a quantum spin model 
with antiferromagnetic (AF) coupling and AF order in one 
direction and ferromagnetic coupling in the perpendicular 
direction with impurities moving through the lattice. The 
impurity motion between sub-lattices couples to quantum 
fluctuations of these pseudo-spin degrees of freedom which 
correspond to the boson hopping from an occupied site of 
the solid to an empty interstitial site. We find that, for 
the limit which describes the case of solid ^i?e, interstitial 
atoms are well-defined delocalized excitations. It is sug- 
gested that during the '^He solid nucleation process, the 
^He impurities stabilize point defects such as dislocations 
or disclinations at the solid-to-liquid interface which be- 
come lines of defects as the solid grows. These defects can 
become highly entangled and can form a glass phase at 
low temperature. The ^He impurities become bound to 
these defects and promote mobile interstitial '^He atoms 
into the disordered solid. We find that such interstitial 
atoms can condense and may be responsible for the ob- 
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served behavior [1]. 

Let us consider a lattice gas model to describe the 
bosonic solid and the added impurities. In such a model, 
we need to consider the interstitial sites as part of the 
lattice and, thus, the ideal quantum solid containing no 
vacancies and no impurities, corresponds to a fractionally 
occupied lattice. For example, the ideal triangular solid 
corresponds to the case of 1/3 filling, namely, to a \/3 x \/3 
ordered solid and the ideal square lattice solid corresponds 
to the \f2 X \f2 checkerboard solid, i.e., 1/2 filling of the 
lattice with bosons. Our model Hamiltonian describing 
a bosonic quantum solid (such as solid ^He) and a small 
concentration of impurities (such as ^He atoms) may be 
written as follows: 

i i 

= ^ [-tbiBlBj + h.c.) + VNiNj] + UY,NI (2) 

<ij> i 

n* =Y, [-Mcj + h.c.) + Whihj] + U*Y, (3) 

<ij> i 

Hint = V*Y1 (^'"J + ^'^^•) + ^' E ^i^i' 
<ij> i 

is the Hamiltonian of the bosonic solid, and Bj (Bi) are 
boson creation (annihilation) operators and Ni = BjBi. 
The fermion operators c|(ci) create (or annihilate) impuri- 
ties on the site i and we have suppressed their spin degree 
of freedom for simplicity. Here, Ui is the impurity number 
operator. We will consider the U ^ oo, U* ^ oo and 
f' — > oo limits (single-site occupation subspace) because 
of the the hard-core interaction. V, V* and W arc positive 
because of the additional energy cost to place an atom in 
an interstitial site. 

In the absence of impurities the well-known pure bosonic 
Hamiltonian (2), in the limit of [/ ^ oo and under 
the transformation B} S+,Bi S~,Ni Sf - \, 

where {Sf , Sf , Sf) arc spin-1/2 operators, reduces to the 
anisotropic spin-1/2 Heisenberg model [12] 

Wb= ^[J5f5; + J.,(5r5j= + 5f5j)]-if^5f, (5) 

<ij> i 

where J = V > and J,j.y = -2tb < and H = fi-z/2V, 
z is the coordination number. 

In order to illustrate the effects of the impurities on 
the stability of the quantum solid and, vice versa, we first 
consider the square lattice because of its simplicity. Fol- 
lowing the general spin-wave (SW) theory for an ordered 
square lattice quantum antifcrromagnct [12], we separate 
the ordered square lattice in two sub-lattices, A (or up, 
or occupied) and B (down, or empty) and we consider bo- 
son operators a] and fej which create spin-deviations with 
respect to the classical Neel ground state in sites of the 
corresponding sub-lattice. The Hamiltonian (5) is approx- 
imated by keeping terms up to quadratic in spin- deviation 



operators; using the Fourier transforms Ok and 6k of the 
operators Oj and bi (as defined in Ref. [12]), where k takes 
values from the Brillouin zone of the \/2 x y/2 sub-lattice 
and introducing the Bogoliubov canonical transformation. 

Ok = UfeCkk + ffe/3^^k' = UkPk - -yfeaLk; where aj^. and 
are boson creation operators, the Hamiltonian takes the 
form 

nE=Eo + Y^ (<4«k + uj^pIp^) , (6) 

k 

where uj^''^ = rfJck ± H, = X'^^l. 7k = 

l/2{cos(h\r) + cosiky)), d = 2 and A = •/.■,,//•/• and 
Uk = v/l/2(l/ek + l), Wk = -s.gn(7k)v/l/2(l/ek- 1). 

It follows from Eq. 6 that for ojq'^ > 0, i.e., for 
H > dJ\/l — A^, the Neel ordered ground state is unsta- 
ble. Ferromagnetic (supcrfluid in the Bosc system) order 
develops in the xy direction and the spins are canted in 
order to acquire a component along the direction of the 
field. The phase diagram obtained for this model with 
this spin-wave approximation agrees reasonably well with 
that obtained by other techniques. [13] When the square 
lattice is half-filled there is a gap G = dJ\^l — A^ for cre- 
ating a propagating pseudo-spin wave excitation, i.e., to 
promote a boson atom to the interstitial band. These in- 
terstitial Quasioarticles move in a band which in the limit 
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Fig. 1: The spectral functions A(k, w) for J/i = 3 (left) and 
J/t = 6 (right) along the Brillouin zone path (0.0) — » (tt, 0) —^ 
(7r,7r)^(0,0). 

We wish to extend this approach to study the motion of 
a single impurity inside the quantum solid. As in Ref. [14], 
we consider as reference state the Neel state with an im- 
purity on the A (substitutional) or B (interstitial) sub- 
lattices. We imagine that there exist operators /ij and 
// which operate on the Neel state and replace respec- 
tively an up-spin or down-spin with an impurity. For 
the single impurity case and by keeping only up to lin- 
ear terms in spin-deviation operators the terms given by 
Eqs. (3) (4) in terms of these operators take the form 

n*+nint = -tj:<ij>,ieAi4f}hi+h.c.)+v*dj:^^sfUi. 



r>.0 



Possible role of ^He impurities in solid '^He 



The entire linearized Hamiltonian (1), using the the Bo- 
goliubov transformation, takes the following form: 



k 



ni 



where ei = 0, £2 = V*d and is given by Eq. (6). 
The operators /j[ and are the Fourier transforms of 

, , (1) , nTTTT. 

-,U. 



fl and h\ respectively and where — 2 /Nt^\^ 



The two Dyson's equations in the 
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Fig. 2: The substitutional impurity moves by means of the 4' ^ 
order process (a)-(e). The interstitial impurity moves via the 
2"'' order process (f)-(h). The impurity, bosonic atoms, and 
empty interstitial sites are denoted as gray solid, black solid, 
and open circles respectively. 
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Fig. 3: The bandwidth of the substitutional and interstitial 
impurity bands as a function of J /t. 



non-crossing approximation are given as follows [14]: 
Gy(k, = , 

- - Eq .9kq G^' (k - q - - CJ^) 



(7) 



where the first equation is obtained for v — 1, v' ~ 2 and 
^ = a and the second equation for j/ = 2, i^' = 1 and 
H = (3. Here a;Q,_^(q) are given after Eq. 6. The Green's 
function Gi (G2) corresponds to the quasi-particles cre- 
ated by hi. ifl). These equations can be solved iteratively 
as in Ref. [14] starting from gI"'' (k, lo) = l/{uj — + irf). 



In our calculations we took V* = V — J , Jxy = 2t 
(i.e., t — tt), H = and e = 0.1. As discussed later 
NMR measurements [15] indicate that for "^iJe impurities 
in solid '^He we should consider J/t >> 1. In Fig. 1 the 
spectral function A(k = 0,uj) is presented for the case of 
J/t = 3 (left) and J/t = 6 (right). Notice that for J/t = 3 
the bandwidth Wa of the substitutional impurity is small. 
For J/t = 6, Wa is very smaU, i.e., WA/t 10~^, and, 
hence, we only show the spectral function of the interstitial 
impurity for this value of J/t (right part of Fig. 1). Notice 
that the spectral function of the interstitial impurity has 
two main peaks, a lower frequency peak with small spec- 
tral weight and a higher frequency one with most of the 
spectral weight. 

In the regime of J >> t, the leading order in t/J which 
allows the substitutional impurity or a ^He atom to move 
is the fourth order process shown in Fig. 2(a-e). Hence, the 
impurity moves in a band with a bandwidth of the order of 
WA/t = A{t/ JY in an expansion of t/J and Jxy/ J (which 
is also small since Jxy = 2i). On the other hand, in our 
case where V* ~V — J , the bandwidth of the interstitial 
impurity is of the order of Wb /t ^ t/ J and it corresponds 
to the process shown in Fig. 2(f-h). States such as those 
of Fig. 2(f) and Fig. 2(h) are connected by second order 
degenerate perturbation theory processes where the states 
shown in Fig. 2(g) are included as intermediate states. 
Namely, the interstitial impurity at site A takes advantage 
of a pair "pseudo-spin" flip near it, as in Fig. 2(g), and 
hops to site B as shown in Fig. 2(h). 

We also carried out a diagonalization in a space which 
includes the 17 states of the type shown in Fig. 2(f-h) and 
their translations through the square lattice. We found 
the dispersion indicated with the black and green lines 
in Fig. 1. The bandwidth of the impurity as calculated 
both by solving the Dyson's Eqs. 7 and by the variational 
approach is shown in Fig. 3 as a function of J/t. Notice 
that the interstitial band (Fig. 1) and bandwidth agree 
very well with those obtained from the Dyson's equations. 
The solid line and dashed fine are fits to Ws/t — B{t/J) 
and WA/t = A{t/JY respectively with A~ B 

Now, let us turn our discussion to the real case of ^He 
impurities in solid '^He. The form for bandwidth of the 
substitutional impurity, i.e., Wa ~ t{t/J)^ (when we take 
i ~ tfe and V* ^ V), is also valid on the triangular and hep 
lattices; namely, when the substitutional band is filled, in 
order for the atoms to move, the same fourth order pro- 
cess, where the atoms momentarily hop to interstitial po- 
sitions, is necessary. In NMR [15] studies tunneling rates 
were found to be of the order of 1 MHz. Using our calcu- 
lated form for the bandwidth Wa — At{t/J)^, and taking 
i = lis: and J = we find Wa ^ x IQ-^K (i.e., 1 
MHz). Using our form for the bandwidth for interstitial 
impurities Wb = Bt^ / J and the same value of t and J 
we find Wb ^ AQmK . Hence, the present theory can re- 
produce (a) the NMR [15] results, (b) that substitutional 
quasiparticles might be localized [11] is solid ^iJe, and (c) 
it suggests that interstitial impuritons might move coher- 
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Fig. 4: Edge-dislocations or disclication linos form as the 
*_ffe solid grows on the cell surface (bottom plane). The top 
(dashed) plane indicates the solid/liquid interface. We imag- 
ine that the solid grows layer by layer. Because the interface 
is two-dimensional, dislocations and disclinations might be ex- 
pected to appear during a relatively fast cooling process at 
any given layer. In addition, the impurities present in the liq- 
uid bind very quickly to such point defects. The point defects 
of the next layer more or less line up with those of the previ- 
ous layer (otherwise it is energetically very costly). Therefore, 
these defects become lines of defects as shown in the above 
figure. When the average distance between such defects is 
smaller than the length Urms which characterizes the wander- 
ing of each defect line, the defect-lines can entangle (such as 
the lines C1C2, D1D2 and E1E2) and can produce a topological 
glass phase. 



ently near the above temperature scale. However, as we 
discuss below, interstitial impurities bind with defects. 

For the case of interstitial ^ile atoms, using our result 
for the bandwidth, W = dJ^y/2J (and J^y ^ 2i), and the 
values of the parameters discussed in the previous para- 
graph, we find W ~ 200mK. This is the temperature 
scale where the possible super-solidity of ^ffe has been 
observed [1]. In the case of real solid ^ife nearest neigh- 
bor empty interstitial sites are much closer to one another 
than occupied sites and, in addition, an interstitial atom 
has higher energy than a substitutional atom and there- 
fore faces a lower potential barrier to tunnel to another 
interstitial site. These two factors increase the tunneling 
frequency by several orders of magnitude relative to the 
observed frequency [15] for substitutional sites. However, 
a finite energy is needed in order to promote atoms to 
the interstitial band [16] and, in addition, vacancies and 
interstitials have a tendency to phase separate in an equi- 
librated crystal [17]. On the other hand, inhomogeneities 
provided by ^He impurities during the ^ife solid nucle- 
ation process or by a relatively fast inhomogeneous cool- 
ing process can create such interstitial atoms as discussed 
next. 

As already discussed, the model studied in this paper is 
very similar to the t — J model used to understand quan- 
tum antiferromagnets, where the impurity degrees of free- 
dom in our case map to the holes in the t — J model and 
the hopping of the bosonic atoms maps to the spin degrees 
of freedom. In the case of the t — J model on the square 
lattice, when a finite density of holes is introduced the 
ground state may have stripe-order [18] where the stripes 



are hole-rich and the domains between stripes are antifer- 
romagnetically ordered (quantum solid domains) with a 
TT-phase shift between domains. As we will argue below, 
a similar phenomenon should be expected to occur in the 
case of an hep lattice of bosons with impurities. Namely, 
edge-dislocations or disclination lines such as those shown 
in Fig. 4 might be created during the process of solid '^He 
growth from the liquid phase. 

First, it is known that the inhomogeneous nucleation of 
solid **ife starts on the cell walls [19]. We will assume 
that solid '^He grows layer-by-layer on the surface of the 
cell for the following reasons: (a) The geometry of the cells 
used in the experiments [1,20] is characterized by a large 
surface with small distances between large confining walls, 
(b) The cooling of the solid is caused through cooling of 
these large surfaces, therefore, atoms near the surface will 
get cold first, (c) Helium atoms at temperature T ~ IK 
get adsorbed on the surface of most substrates because of a 
relatively strong dipolar interaction with the surface [21]. 
In studies oi'^He on graphite [22,23] for example and in all 
known substrates the first layer of adsorbed *_ffe is solid. 
The second layer of '^He on the graphite surface is also 
solid at layer completion density [22,24]. After the first 
or a few layers are deposited on the surface, the adsorbed 
solid '^He layer is compressed relative to the liquid in con- 
tact and, therefore, the coated surface provides a stronger 
attractive potential for an atom as compared to that pro- 
vided by its surrounding atoms in the liquid; thus, the 
atoms tend to get adsorbed and form a new solid layer. 
Nucleation of solid '^He had been a puzzle in the past just 
because of the large surface energy cost required in homo- 
geneous nucleation [19]. In a layer-by-layer growth this 
is not an issue, because when the next layer is deposited, 
there is no additional surface area introduced, namely, just 
the same-size interface advances. 

It is believed [25] that the solid/liquid interface in the 
presence of ^He impurities becomes rough due to bind- 
ing of such impurities to the interface. Even an extremely 
low impurity concentration, namely as low as 1 ppm, has 
a significant effect in this process [25]. In this paper, it 
is proposed that the ^He impurities bind to defects [26], 
namely to dislocations or disclinations, which are expected 
to form on the two-dimensional solid/liquid interface [27]. 
There are various theories of two dimensional melting [28], 
the most popular of which is the theory of Halperin, Nel- 
son and Young [27] which was inspired by the Kosterlitz- 
Thouless(KT) theory of vortices in superfluid films. Ac- 
cording to this theory, the melting proceeds via two contin- 
uous KT transitions, the lowest temperature one caused 
by the unbinding of dislocations and a higher tempera- 
ture transition caused by the unbinding of disclinations. 
Both types of defects are point-like singularities; the for- 
mer corresponds to the singular behavior of the atomic 
lattice displacement field and the latter is associated with 
the singular behavior of the angle field 9{r) which char- 
acterizes the fluctuations in the bond orientational order. 
Our following analysis applies equally to both type of de- 
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fects and, therefore, we will use the term defect in general. 
As the solid/liquid interface is cooled down at a relatively 
fast rate, locally small-size 2D clusters of atoms at the 
solid/liquid interface achieve local crystalline order rather 
quickly; however, the time scale to achieve a global equi- 
librium is much longer because these clusters may have 
different orientations. In order to describe the state of 
this globally disordered but locally ordered 2D system, 
we need to imagine the presence of dislocations (or discli- 
nations). Generally, as a 2D disordered solid is cooled 
down to low temperature, pairs of such thermally created 
excitations having opposite topological charge move very 
slowly towards each other and are annihilated (or bind). 
However, the ^He impurities which are present in the liq- 
uid find it energetically favorable to bind to these defects 
and the combined impurity/defect system becomes stable. 
Furthermore, it becomes energetically favorable for these 
lattice defects of the next deposited solid layer to follow 
the positions of the defects on the previous layer, thus, 
defect-lines are formed as the interface advances; namely, 
lines of singularities grow perpendicular to the interface 
(see Fig. 4). 

The root-mean-square projection Ums of the end-to-end 

vector of a defect-line is given as [29] u^ms — V < > — 
{2TrkBTd/ey^^ where d is the thickness of the solid film 
and e is the line tension along the defect-line. When Urms 
is a few times greater than the average distance between 
such defects, neighboring defect-lines will entangle. Tak- 
ing T ~ liiT and e ~ lOK/A and d = 0.3mm we find 
that iirms ^ 10^ A which is greater than the average dis- 
tance between impurities and, thus, entanglement seems 
likely. When the solid is cooled to low temperature defects 
with opposite topological charge cannot be mutually an- 
nihilated because the entangled defects cannot move. The 
entanglement is topologically protected and this could cre- 
ate a topological glass state of solid ^iJe. This is analo- 
gous to the vortex-glass state proposed for superconduc- 
tors [30,31]. 

As discussed earlier '^He impurities, even at concentra- 
tions as low as 1 ppm, infiuence very significantly the solid 
^He nucleation process [25]. In particular even such very 
small amount of impurities can change the roughening 
transition temperature [25] by 20%. This has been in- 
terpreted as the result of ^He impurity adsorption at the 
liquid to solid interface. As discussed such impurities sta- 
bilize the defects and as the interface advances defect lines 
form which may entangle and this can influence signifi- 
cantly the roughening transition. Therefore, it is possible 
that the ^He impurities, through the creation of such de- 
fects, promote ^He atoms to the interstitial band which 
can move from domain to domain and by means of the 
entangled defects. This creates a metastable disordered 
supersolid in which the carriers flow through a topolog- 
ical glass of entangled defect-lines. This scenario might 
explain the observed hysteresis [2] by annealing of the 
solid helium sample. Most recently Rittner and Reppy [20] 
have observed very high superfluid response in their tor- 



sional oscillator experiments of solid '^He samples grown 
from the liquid phase by very rapid cooling. Furthermore, 
these samples have high surface-to-volume-ratio and high 
concentration of such defects is expected to occur under 
such conditions. Namely, under rapid cooling such defects 
are expected to occur in the 2D solid/liquid interface, be- 
cause, the atoms order locally to form a microscopic-size 
solid very quickly, but the time scale for the annihilation of 
pair of defects of opposite topological charge is very long, 
namely much longer than the time required to establish 
local (i.e., at a microscopic scale) equilibrium. Once these 
defect-lines are formed, they wander around as shown in 
Fig. 4 and they may entangle around each other several 
times because of their high density. The process of an- 
nihilation of these defects is then very difficidt because 
of their entanglement and this may lead to an extremely 
long-lived topological glass phase. 

Experimentally, while the critical temperature in- 
creases with increasing ■^He impurity concentration X3, for 
the superfluid response ps there is an optimum X3 above 
which Ps decreases with increasing X3. This behavior is 
consistent with the theory presented here. The degree 
of inhomogeneity is proportional to the density of defects 
which is also proportional to the impurity concentration 
Xs. The carrier density, namely the density of intersti- 
tial atoms, increases with X3, therefore, is expected to 
increase with X3. Since Tc is more or less a locally de- 
termined quantity, the efi^ect of global disorder on Tc is 
much less important than the fact that the number of car- 
riers rise with disorder. The superfluid response, however, 
is expected to depend strongly on phase fluctuations of 
the superfluid order parameter, therefore, disorder is ex- 
pected to have harmful effects to the long-range coherence. 
In summary, it is reasonable to expect that at high con- 
centration of such ^He impurities, the disorder which is 
caused by the lattice defects should significantly harm the 
long-range coherence, however, Tc may continue to rise be- 
cause, while the long-range disorder does not significantly 
influence the value of Tc, Tc increases due to the increase 
of the number of carriers. 

I would like to thank K. Wicrschem for proof-reading 
the manuscript and M. Boninsegni for useful discussions. 
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Possible role of '^He impurities in solid '^He 
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We use a quantum lattice gas model to describe the essential aspects of the motion '^He atoms 
and of a ^He impurity in solid *//e. We find that '^He impurities promote ^He atoms to interstitial 
sites and this can turn the bosonic quantum crystal into a metastable supersolid. It is suggested 
that ^He impurity atoms, which produce the interstitial *ffe atoms, might have been reabsorbed 
by pure "^He solid formed during the first stage of a multi-stage solid ^He nucleation process. While 
we use the "spin"-wave approximation and low dimensional lattices to illustrate some of the ideas, 
we argue that the conclusions drawn from these studies may be valid for the real system. 



Kim and Chan[il,[i] (KC) usine; the torsional oscillator 
technique found a decrease in the resonant period of solid 
He confined in porous vycor glass Q and in bulk solid 
helium below 2QQmK, indicating the possible onset of 
superfluidity in solid helium. The experimental results of 
KC have been independently confirmed by Rittner and 
ReppyQ (RR) using a different geometry to confine the 
solid. In addition, RR observed that the superfluidity of 
solid helium can be significantly influenced or eliminated 
by annealing of the solid helium sample. Because of these 
history-dependent results and the negative results in at- 
tempts to drive flow by pressure^], the interpretation of 
the results of KC is subject to debate. 

Furthermore, it has been observed that the superfluid 
response and superfluid fraction are strongly dependent 
on the amount of isotopic ^He impurities which exist in 
the naturally available helium. When, for example, the 
naturally occurring concentration of ^He impurities of 
about one partper million is reduced to less than one 
part per billion the superfluid fraction is reduced from 
about 1% to approximately 0.03%. When the ^He impu- 
rity concentration is increased to about 0.1% the super- 
fluid fraction vanishes^. Previously, Ho et al^ found 
an anomalous behavior of the acoustic attenuation and 
velocity in solid '^He below ^ 2QQmK at a low concen- 
tration of ^He impurities. 

The possibility of superfluidity of solid '^He has been 
extensively discussed 0,11, @, 0, Eil . It is of fundamental 
value for our understanding of this and related phenom- 
ena because of the implication that there can be coexis- 
tence between spatial and momentum space order [l2j. 

In this paper we study the role of very low density of 
impurities in a bosonic hard core solid. While the the- 
oretical studies of the role of impurities in solid helium 
has a long history 0, [H, 14 1, the present study, based on 
an analogy with models of doped quantum antiferromag- 
nets and by using rather recently developed techniques 
for such systems, sheds some "new" light on the problem. 
Our model, that describes the impurity motion in an oth- 
erwise ideal quantum bosonic crystal, maps to a quantum 
spin model with antiferromagnetic (AF) coupling and AF 
order in one direction and ferromagnetic coupling in the 



perpendicular direction with impurities moving through 
the lattice. The impurity motion between sub-lattices 
couples to quantum fluctuations of these pseudo-spin de- 
grees of freedom which correspond to the boson hopping 
from an occupied site of the solid to an empty interstitial 
site. We find that substitutional impurities become lo- 
calized and an interstitial impurity relaxes into a substi- 
tutional site by promoting a '^He atom to the interstitial 
band. Furthermore, substitutional ^He impurities, be- 
cause of their larger zero point motion, tend to promote 
^He atoms into the interstitial band. We discuss that the 
bandwidth of the interstitial band is of the same order 
of magnitude as the temperature scale below which the 
non-classical moment of inertia was observed [J, 

Let us consider a lattice gas model to describe the 
bosonic solid and the added impurities. In such a model, 
we need to consider the interstitial sites as part of the 
lattice and, thus, the ideal quantum solid containing no 
vacancies and no impurities, corresponds to a fraction- 
ally occupied lattice. For example, the ideal triangular 
solid corresponds to the case of 1/3 filling, namely, to a 
\/3 X -\/3 ordered solid and the ideal square lattice solid 
corresponds to the x y/2 checkerboard solid, i.e., 1/2 
filling of the lattice with bosons. Our model Hamilto- 
nian describing a bosonic quantum solid (such as solid 
^He) and a small concentration of impurities (such as 
^He atoms) may be written as follows: 

7i: = 7iB + ?i:*+7i:,„t-Ai^iv,-/i*$^"^, (i) 

i i 

Ub^Y. i-tbiBjB, + h.c.) + VN,N,] + uY,Nf, (2) 

<ij> i 

n* ^Y. [-*(4cj + h.c.) + Wh,hj] + U*Y,f^l (3) 

<ij> i 
<ij> i 

TLb is the Hamiltonian of the bosonic solid, and B\{Bi) 
are boson creation (annihilation) operators and Ni — 
BjBi. The fermion operators c|(ci) create (or annihilate) 
impurities on the site i and we have suppressed their spin 



2 



degree of freedom for simplicity. Here, rii is the impurity 
number operator. We will consider the J7 — > oo, ?7* ^ oo 
and U' ^ CO Hmits (single-site occupation subspace) be- 
cause of the the hard-core interaction. V, V* and W are 
positive because of the additional energy cost to place an 
atom in an interstitial site. 

In the absence of impurities the well-known pure 
bosonic Hamiltonian ([2]), in the limit of U —^ oo and un- 
der the transformation Bj , Bi S^^ ,Ni Sf — ^, 
where {Sf,Sf, Sf) are spin-1/2 operators^ reduces to the 
anisotropic spin-1/2 Hcisenberg model[15| 

Tis = ^ [JStS^^ + J^yiS^S^ + 5f5j)] -HY^St, (5) 

<i3> i 

where J = V > Q and J^y = —2th < and H = fi—z/2V, 
z is the lattice site coordination number. 

In order to illustrate the effects of the impurities on 
the stability of the quantum solid and, vice versa, we 
first consider the square lattice because of its simplic- 
ity. Following the general spin-wave (SW) theory for an 
ordered square lattice quantum antiferromagnetTs'l, we 
separate the ordered square lattice in two sub-lattices, 
A (or up, or occupied for the case of the solid) and B 
(down, or empty) and we consider boson operators a| 
and &| which create spin-deviations with respect to the 
classical Neel ground state in sites of the correspond- 
ing sub-lattice. The Hamiltonian ^ is approximated 
by keeping terms up to quadratic in spin-deviation op- 
erators; using the Fourier transforms ay and 6k of the 
operators Ui and bi (as defined in Ref. llSh , where k takes 
values from the Brillouin zone of the V2 x \/2 sub-lattice 
and introducing the Bogoliubov canonical transforma- 
tion, Ok = UfcOfk + VkP'Li^, = UkPk - I'fca^k' where 
Qfj^ and are boson creation operators, the quadratic 
linearized Hamiltonian can be diagonalized and it takes 
the form 



Eo + J2{uj^ala^ + u;^Pll3^y (6) 

k 



where lu^'^ = dJe^ ± Ck = ^/l 
l/2{cos{kx) + cos{ky)), d = 2 and A 



A^T^, 7k = 
Jxy/J, pro- 
vided that we choose the functions Wk and Wk as follows: 
Mk = v/l/2(l/ek + l), Vk = -S5"(7k)VV2(l/ek- 1). 

It follows from Eq. [6] that for uJq'^ > 0, i.e., for 
H > dJ\/l — A^, the Neel ordered ground state is unsta- 
ble. Ferromagnetic (superfluid in the bose system) order 
develops in the xy direction and the spins are canted in 
order to acquire a component along the direction of the 
field. The phase diagram obtained for this model with 
this spin-wave approximation agrees reasonably well with 
that obtained by other techniques. [l6| When the square 
lattice is half-filled there is a gap G = dJy/1 — A^ for 
creating a propagating pseudo-spin wave excitation, i.e., 
to promote a boson atom to the interstitial band. These 



interstitial quasiparticles move in a band which in the 
limit of Jxy « J has a bandwidth W = dJ^y/2J. 

We wish to extend this approach to study the motion of 
a single impurity inside the quantum solid. As in Ref. 1 171. 
we consider as reference state the Neel state with an im- 
purity on the A (substitutional) or B (interstitial) sub- 
lattices. We imagine that there exist operators /i| and // 
which operate on the Neel state and replace respectively 
an up-spin or down-spin with an impurity. For the single 
impurity case and by keeping only up to linear terms in 
spin-deviation operators the terms given by Eqs. p]) (|4]) 
in terms of these operators take the form H* + Hint — 

-iE<.,>,eA(«l/>. + h.c.) + v*dY.,^^flu 

The entire Hamiltonian ([1]) when linearized and for the 
case of a single impurity, as described above, and by using 
the Bogoliubov transformation for the boson operators, 
takes the following form: 



k,q 



n 



L : 



where ti = 0, 62 = V*d and is given by Eq. 
The operators and are the Fourier transforms of 



fl and h\ respectively and where = d\/2/N tjy^ 



9uo = d,/2/Nt^^v. 



The two Dyson's equations in 




FIG. 1: The spectral functions A(k,w) for J/t = 3 (left) and 
J/t — 6 (right) along the Brillouin zone path (0.0) (n, 0) 
{7v,t:) ^ (0,0). 



the non-crossing approximation are given as follows 17| : 

1 



- - Eq 5kq ^Gi,' (k-q-uj-uj^)' 



(7) 



where the first equation is obtained for = 1 , i^' = 2 and 
^ — a and the second equation for j/ = 2, z/' = 1 and fi = 
(3. Here uja,0{q) are given after Eq. [S] The Green's func- 
tion Gi (G2) corresponds to the quasi-particles created 
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FIG. 2: The substitutional impurity moves by means of the 
4"* order process (a)-(e). The interstitial impurity moves via 
the 2""* order process (f)-(h). The impurity, bosonic atoms, 
and empty interstitial sites are denoted as gray solid, black 
solid, and onen circles resoectivelv. 
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FIG. 3: The bandwidth of the substitutional and interstitial 
impurity bands as a function of J /t. 



by if I)- These equations can be solved iteratively as 
in Ref . [13 starting from ci"^ (k, u) — 1 / {uj — e^j + irf) . 

In our calculations we took V* ^ V = J, Jxy = 2t 
(i.e., t = tb), H = and e = 0.1. As discussed later 
NMR measurements [3| indicate that for '^He impurities 
in solid '^He we should consider J/t >> 1. In Fig. [T] 
the spectral function A(k — 0, uj) is presented for the 
case of J/t = 3 (left) and J/t — 6 (right). Notice that 
for J/t = 3 the bandwidth Wa of the substitutional im- 
purity is small. For J/t = 6, Wa is very small, i.e., 




FIG. 4: 



A substitutional '^He impurity in solid *i?e pro- 



motes He atoms from the substitutional to the interstitial 
band. The size of the black circles is proportional to the *He 
occupancy probability. The white circles indicate empty in- 
terstitial sites. The site A is occupied by a '^He impurity. 



WA/t ~ 10~^, and, hence, we only show the spectral 
function of the interstitial impurity for this value of J/t 
(right part of Fig. [T]). Notice that the spectral function 
of the interstitial impurity has two main peaks, a lower 
frequency peak with small spectral weight and a higher 
frequency one with most of the spectral weight. 

In the regime of J >> t, the leading order hit/ J which 
allows the substitutional impurity or a '^He atom to move 
is the fourth order process shown in Fig. [IJa-e). Hence, 
the impurity moves in a band with a bandwidth of the or- 
der of WA/t = A{t/ JY in an expansion oft/ J and Jxy/ J 
(which is also small since J^y = 2t). On the other hand, 
in our case where V* — V ^ J, the bandwidth of the 
interstitial impurity is of the order of Ws/t ~ t/J and 
it corresponds to the process shown in Fig. El^f-h). The 
states such as those of Fig.[2ljf) and those of Fig. [2th) are 
degenerate and are connected by second order degenerate 
perturbation theory processes where the states shown in 
Fig. Wis) are included as intermediate states. Namely, 
the interstitial impurity at site A takes advantage of a 
pair "pseudo-spin" flip near it, as in Fig. [S^g), and hops 
to site B as shown in Fig. [DJh). 

We also carried out a diagonalization in a space which 
includes the 17 states of the type shown in Fig.[2l^f-h) and 
their translations through the square lattice. We found 
the dispersion indicated with the black and green lines 
in Fig. [TJ The bandwidth of the impurity as calculated 
both by solving the Dyson's Eqs.[7]and by the variational 
approach is shown in Fig. [3] as a function of J/t. Notice 
that the interstitial band (Fig. [1]) and bandwidth agree 
very well with those obtained from the Dyson's equations. 
The solid line and dashed line are fits to Ws/t = B{t/J) 
and WA/t A{t/J)^ respectively with ^ ~ S ~ 1. 

Now, let us turn our discussion to the real case of 
^He impurities in solid '^He. The form for bandwidth 
of the substitutional impurity, i.e., Wa ^ t{t/J)^ (when 
we take t ^ t^ and V* ^ V), is also valid on the trian- 
gular and hep lattices; namely, when the substitutional 
band is filled, in order for the atoms to move, the same 
fourth order process, where the atoms momentarily hop 
to interstitial positions, is necessary. In NMR[i8| stud- 
ies tunneling rates were found to be of the order of 1 
MHz. Using our calculated form for the bandwidth 
Wa = At{t/Jf, and taking t = IK anA J ^ 30if, we 
find Wa 4 x IQ'^K (i.e., 1 MHz). Using our form for 
the bandwidth for interstitial impurities Wb = Bt^/J 
and the same value of t and J we find Wb ~ AOmK. 
Hence, the present theory can reproduce (a) the NMR [lit 
results, (b) that substitutional quasiparticles might be 
localized|14| is solid '^He, and (c) it suggests that inter- 
stitial impuritons might move coherently near the above 
temperature scale. However, as we discuss below inter- 
stitial impurities lighter than the '^He atoms very quickly 
relax into substitutional sites. 

For the case of interstitial ^77e atoms, using our result 
for the bandwidth, W = dJ^y/2J (and J^y ~ 2t), and the 
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values of the parameters discussed in the previous para- 
graph, we find W ~ 200mK. This is the temperature 
scale where the possible super-solidity of ^He has been 
observed [HQ- In the case of real solid *He nearest neigh- 
bor empty interstitial sites are much closer to one another 
than occupied sites and in addition, an interstitial atom 
has higher energy than a substitutional atom and there- 
fore faces a lower potential barrier to tunnel to another 
interstitial site. These two factors increase the tunneling 
frequency by several orders of magnitude relative to the 
observed frequency[l8| for substitutional sites. However, 
an energy gap is needed in order to promote them to the 
interstitial band and the role of '^He impurities becomes 
important right here as discussed next. 

An interstitial impurity is expected to quickly relax to 
a substitutional site by promoting a '^He atom to an in- 
terstitial position. The reason is that the ■^He atomic 
mass is smaller than the '^He mass, therefore, the con- 
figuration in which the ^He atoms is "squeezed" in the 
interstitial space has significantly higher energy than the 
one where the ^iJe is in the interstitial space. For a simi- 
lar reason even a substitutional ^He atom has a tendency 
to promote a "^He atom to an interstitial space as illus- 
trated in Fig. m A substitutional ^He impurity in solid 
^iJe prefers a larger room for zero-point motion than the 
solid "^He lattice spacing because of its lighter mass. As 
a result, it causes a local expansion which means that 
the atoms labeled as 1-6 in the Fig. [4] will spend some 
time as interstitial '^He atoms. This implies that ^He 
impurities are expected to promote '^He atoms from the 
substitutional band which is full to an empty interstitial 
band with a relatively large bandwidth. 

Interstitial impurities and interstitial "^He atoms could 
be created by the non-equilibrium nucleation process that 
produces solid *_ffe from the liquid phase. The kinetics 
of nucleation and of phase separation of solid ^He He 
mixtures is still under experimental investigation 3, 23, 
21| . It is known that there is a great asymmetry between 
the kinetics of phase separation and homogenization in 
solid ^He — ^ He mixtures. The observation of anoma- 
lously fast mass transfer at inclusion-dissolution of ^He 
in solid ^He matrix was explained by ballistic transfer 
of ^He quasiparticles[2lj. "^iJe requires a lower pressure 
than ^He to solidify; thus, it is possible that, at a first 
stage, the nucleation of pure '^He takes place and the ^He 
impurity atoms are excluded when solid clusters form. At 
a later stage these ^He atoms may be slowly reabsorbed 
through diffusion inside the solid '^He crystal, which, as 
discussed, produce interstitial ^iJe atoms which should 
remain as metastable defects in the crystal. This scenario 
might explain the observed hysteresis Q by annealing of 
the solid helium sample. 

If indeed '^He impurities drive the superfluid behavior 
in solid ^He, the condensate fraction, naively, should be 
proportional to the '^He impurity concentration. There- 



fore, the large difference between condensate fraction and 
the observed superfluid fraction needs to be explained. A 
large difference between these cannot be excluded be- 
cause we have examples, such as liquid '^He near so- 
lidification and at low temperature, where the conden- 
sate fraction is small and the superfluid fraction is al- 
most 100%. However, a simple explanation from the 
present work cannot be provided and further experimen- 
tal and theoretical work is needed. Some of the pro- 
posed ideas can be tested using path integral Monte Carlo 
simulationj22|. 
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